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ABSTRACT 


The  equations  of  motion  of  a  beam  with  the  lateral 
surface  thermally  insulated  are  derived,  including 
the  effects  of  shear  deformation  and  rotatory  iner¬ 
tia.  Thermoelastic  coupling  in  the  heat  conduction 
equation  and  in  the  elastic  constitutive  relations  is 
also  included.  Axial  deformations  of  the  beam  are 
taken  Into  account.  Two  examples  are  presented 
and  discussed. 
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SYMBOLS 


A 

C 


F(t),  G(t) 

I(w) 

^xx*  lyy’  *xy 
Kx*  Ky 


Mx.  My 
N 


Px*  Py’  Pz 

Qx'  Qy 

Rx'  Ry 
T 

T 

To 

Tx*  Ty 

U.  V.  w,  Wx,  wy 

u,  v,  w 

•a,  v,  w 


cross-sectional  area 

boundary  of  cross  section  of  a  beam 

specific  heat  at  constant  volume 

E/p 

modulus  of  elasticity 

u  +  v  +  W  „ 

,  x  ,  y  ,  z 

see  pages  18  and  20 
imaginary  part  of  w 
moments  of  inertia 
Timoshenko  shear  constants 
thermal  conductivity 
bending  moments  (see  page  5) 
axial  force  (see  page  5) 
mode  number 

applied  forces  per  unit  length  (see  page  5) 

shear  resultants  (see  page  5) 

applied  moments  per  unit  length  (see  page  5) 

/ATdA 

temperature  measured  from  ambient 
ambient  temperature  (absolute) 
fA  y,  x  TdA 

displacement  resultants  (see  page  5) 

l/A^u.v, w  dA 

displacements 

iv 


x,  y.  z 


a 

V 


r. 


*!•  **2’  ^ 3 

V 

l 

K 


c  ,  <r  ,  etc. 
xx  yy 


0  »  0 
rx  y 


rectangular  Cartesian  coordinates 
coefficient  of  thermal  expansion 
fcj/(«nir)  (see  page  19)> 

|n2ir2/(*f2)[EIyy/(pA)]  1  ^ 2 J/ ( A/Iyy  +  n2Tr2/f  2)(see  page  22) 

|(3\  +  2p)a2T0]/[(X  +  p)pc] 
wl  / (nirc  j )  (see  page  19), 

wl2/(n2ir2)[pA/(EIyy)]1/2  (see  page  22) 

roots  of  the  frequency  equations,  Eqs.  (39)  and  (49) 

Poisson's  ratio 

<-i>l/2 

k/pc  thermal  diffusivity 
Lame  elastic  constants 
density 

stress  components 

cross-sectional  rotations  (see  page  13) 
frequency 
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SECTION  I 


INTRODUCTION 

In  treating  the  problem  of  calculating  atreaaea  and  deflections  of  a 
complicated  structure,  the  common  practice  is  to  consider  the  structure  as 
being  made  of  a  number  of  simpler  structural  elements  whenever  possible 
and  to  treat  these  elements  separately.  Of  all  of  the  simple  structural  ele¬ 
ments,  the  beam  is  perhaps  the  most  common  and  at  the  same  time  the  most 
useful  element.  In  many  applications,  the  beam  can  be  thought  of  as  having 
dynamic  tractions  and  temperature  gradients  applied  to  its  bounding  surfaces. 
This  report  treats  the  problem  of  the  dynamic  response  of  a  beam  to  dynamic 
loads  and  to  thermal  gradients  applied  to  the  ends  of  the  beam  under  the  con¬ 
dition  that  the  lateral  surface  of  the  beam  is  insulated.  This  latter  condition 
could  be  fulfilled  approximately  if  the  beam  were  in  the  hard  vacuum  of  a 
space  environment  and  if  the  thermal  radiation  from  its  lateral  surfaces  could 
be  neglected. 

1  2 

It  is  well  known  ’  that  the  temperature  couples  with  the  elastic  prob¬ 
lem  in  two  ways:  the  coupling  in  the  constitutive  relations  described  by 

3  1 

Neumann  and  the  coupling  in  the  heat  conduction  equation  described  by  Biot. 

2 

This  latter  coupling  can  be  neglected  for  a  wide  class  of  problems,  and 

several  articles  have  treated  the  problems  associated  with  this  assumption.*  ^ 

7 

Ignaczak  and  Nowacki  have  presented  a  similar  analysis;  however,  because 
the  reference  was  not  available  to  the  author,  no  appraisal  of  it  relative  to 
the  present  work  could  be  made.  The  Biot  coupling  is  included  in  this  report 
in  a  general  theory  which  also  includes  rotatory  inertia  and  shear  deformation. 
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thermal  conductivity;  p  is  the  density;  T  is  the  temperature;  c  is  the  specific 

heat  of  constant  volume;  e  is  the  dilatation  and  is  given  by  u  +  v  +  w 

» y  #  * 

A  comma  followed  by  a  subscript  denotes  differentiation. 

The  terms  (3\  +  2p)o  T  occurring  in  Eq.  (3)  are  the  usual  Neumann 
i  8 

coupling  terms.  ’  The  term  ( 3 X.  +  2p)a  Tn  e  in  Eq.  (3)  is  the  coupling 
j  u  ,  t 

described  by  Biot.  The  method  of  procedure  is  to  integrate  Eqs.  (1)  through 
(4)  across  the  thickness  of  the  beam.  The  centroid  of  the  cross  section  of  the 
beam  lies  along  the  z  axis,  as  shown  in  Fig.  1.  Green's  lemma  is  used  in 
the  form 


+  Q  JdA 


=  /• 


Q  dx  +  P  dy 


where  A  denotes  the  cross-sectional  area  of  the  beam  and  C  the  boundary  of 
the  cross  section.  The  equations  will  be  derived  first  neglecting  the  shear 
deformation  of  the  beam  which  will  be  put  in  later. 
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The  phyaical  meaning  of  these  quantities  is  clear:  Q  and  Q  are  the  shear 

x  y 

forces  and  N  is  the  tensile  force  in  the  beam;  Mx  and  are  the  usual  bending 

moments;  U,  V,  W,  W  ,  W  are  the  displacement  resultants  and  moments 

x  y 

and  will  be  discussed  later;  px,  p^,  p2>  are  the  applied  pressures,  taken 
positive  when  acting  in  the  direction  opposite  to  the  respective  coordinate 
areas;  R  and  R  are  the  applied  bending  moments,  taken  positive  in  the 

x  y 

aforementioned  sense. 

Equations  (1)  are  now  integrated  over  the  cross  section  of  the  beam. 

The  last  of  Eqs.  (1)  is  multiplied  by  x  and  y  respectively  and  integrated.  The 
results  are: 


Qx.z- 

Px  = 

pu,tt 

Qy. z  " 

py  = 

pv.tt 

N  * 

»  36 

pz  = 

pw,tt 

M  -  Q„  - 

y.  36  X 

R  = 

y 

pWy.  « 

M  -  Q,  - 

X,  z  y 

Rx  = 

PWx,  tt 

The  same  procedure  is  now  applied  to  the  heat  conduction  equation, 

Eq.  (2).  Integrating  Eq.  (2)  across  the  thickness,  one  has 

“  fA  <T.  XX  +  * .  y„'d A  zzdA-  <"=A.  ,dA  * 

(3\  +  2p)a  Tq  /  e  tdA  (8) 
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If  Green's  lemma  is  used  and  if  it  is  assumed  that  differentiations  with  respect 
to  z  and  t  can  be  interchanged  with  the  integration,  one  has 

k  T, -  P«  T,  t  +  kL{-  T.  y*1  +  T.  *dy>  ■ 


where 


(3\  +  2^)a  T Q£f  e  dA 


(9) 


T  =  /  T  dA 
"A 


(10) 


If  the  lateral  surface  of  the  beam  is  assumed  to  be  insulated,  then  the  integral 

( -  T  dx  +  T  dy)  vanishes  since  the  integrand  is  proportional  to  8T/8n. 

C  ,  y  >  x 

Then  (9)  becomes 


/c 


kT.z,-"c  T.  .  5  <3X  +  T0 


(11) 


Multiplying  Eq.  (2)  by  x  and  y  respectively  and  integrating  yields 


k  Ty,  sz  *  Pc  Ty,  t  *  kfA  T,  xdA  =  (3K  +  2^)a  T0  It  fA  Xe  dA 

k  Tx.  zz  ‘  Pc  Tx,  t  *  k  /  T,  ydA  =  (3X  +  T0  H  L  yC  dA  J 


(12) 


(  ) 
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where 


T 


x 


(13) 


This  is  as  far  as  one  can  go  without  assuming  some  sort  of  deformation 
and  temperature  variation  across  the  thickness  of  the  beam.  For  the  defor- 
mat  ions,  it  is  not  unreasonable  to  expect  a  linear  variation;  this  will  be  taken 
to  be 


u  =  u(z,  t)  ,  v  =  v(z,  t) 

> 

w  =  w(z,  t)  -  xu  -  yv 

9  z  ,z 


(14) 


This  assumption  corresponds  to  Bernoulli- Euler  bending  theory.  To  include 
the  effects  of  shear  deformation,  this  assumption  must  be  modified  slightly; 
this  point  will  be  taken  up  in  detail  later.  Since  the  outer  surface  of  the  beam 
is  insulated,  it  would  not  be  unreasonable  to  assume  a  linear  variation  in 
temperature  as  well.  For  other  boundary  conditions,  such  as  a  given  tempera¬ 
ture  on  the  surface,  this  assumption  would  not  necessarily  be  warranted. 

With  the  assumption  of  linear  temperature  gradients,  and  using  the 
definitions  (10)  and  (13),  one  can  express  the  temperature  as 


T=  Tiz^ti  + 


I  T  (z,  t)  -  I  T  (z,t) 
xx  y  xy  x 

I  I  ^~r 

xx  yy  xy 


+  y 


I  T  (z,  t)  -  I  T  (z,t) 
yy  x'  1  xy  y*  ' 

i  i  -  r 
xx  yy  xy 


(15) 
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f 


o 


?  ( 


i 


1 


This  assumption  gives  rise  to  a  fundamental  inconsistency.  The  boundary 
condition  for  an  insulated  surface  states  that  the  normal  derivative  of  the 
temperature  vanishes  on  the  lateral  surface  of  the  beam.  Equation  (15)  will 
not,  in  general,  satisfy  this  condition.  It  is  felt,  however,  that  this  incon¬ 
sistency  is  an  approximation  of  accuracy  comparable  to  the  usual  approxima¬ 
tions  of  beam  theory,  and  should  yield  a  good  degree  of  quantitative  accuracy. 
A  similar  assumption  was  made  by  Biot*  in  a  similar  problem  which  was 
solved  by  a  variational  procedure. 

The  assumptions  (14)  imply  that 


<r  =<r  =0 

xx  yy 


<r  =  «r  =<r_  =  0 

xy  yz  xz 


ffzz  =  E(w,  z  •  xu,  zz-yv,zz>‘EaT 


aT 


e  =  (1  -  2v)(wi  z  -  xu  zz  -  yv  M)  +  (3\  +  2p)  ^ 


(16) 


These  implications  are  inconsistent  with  the  assumptions  used  in  deriving 
Eqs.  (7).  These  inconsistencies  are  inherent  in  beam  theory,  and  the  justi¬ 
fication  of  the  deformations  as  given  by  Eqs.  (14)  has  been  established  by 
centuries  of  use. 

The  stress  resultants  become 


o 


N  =  EAw  _  -  EaT 

EaTy 

EaTx 
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M  =  -E(I  u  +  I  v  J 
y  '  yy  ,zz  xy  ,  zz' 


M  =  -E(I  u  +  I  v  ) 
x  '  xy  ,  zz  xx  ,zzr 


(17) 


From  the  third  equation  of  (16),  it  is  easily  shown  that 


KI  IM-IM  I  M  -  I  M 
N  +  x  xx  y  V  x  +yJQL_5 - *L-J 1 

zz  *  ~ i  ii  -i2 

xx  yy  xy  xx  yy  xy 


I  T  -  I  T  I  T  -  I  T 

Eox  J2UL—*y  *  .  Eay  JT-J?  *Y  Y  .  Eal 


II  -  I2 

XX  yy  xy 


(18) 


I  I 
xx  yy  xy 


Using  the  last  of  Eqs.  (16)  plus  Eq.  (15)  to  calculate  8T/8y  and  8T/8x,  one 
obtains 


kT,zz*Pc 


1  - 


(3\  +  2*i)Vt0 


t 


(3\  +  2p)a  Tq(1  -  2v)Aw>zt 


k  T  -  pc 
y.  zz  K 


1 


(3\  +  2p)a*T0 


I  T  -  I  T 

T  t  -  k  A  xx  = 

y*  I  I  -  r 

xx  yy  xy 


(\  +  p)pc 

-(3\  +  2p)o  Tn(l  -  2v)(I  u  *  +  I  v  .) 
'  r'  0  yy  ,  zzt  xy  ,  zzt' 


k  Tx,  zz  -  Pc 


(3\  +  2p)2Q2T0 


(X.  +  p)pc 


I  T  -  I  T 
T  -  kA-BL_* - = 


y.t 


II  -  I 
xx  yy  xy 


.(3X  +  2p)aT0(l-2v)(Ixyu  zat  +  I^v 


(19) 
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It  now  remain*  only  to  apply  assumptions  (14)  to  the  equations  of 
motion,  Eqs.  (7).  The  result  is 


Qx,Z-Px  =  PAu,  tt 


Qy,  z  *  Py  *  PAv,  tt 


N,  z-pz  =  pAw  u 


^y,  z  •  Qx  '  Ry  p^yyU>  2tt  +  *xyv>  ztt^ 


MX,  z  -  Qy  -  Rx  -  “P^xyu,  ztt  +  *xxV,  ztt^ 


(20) 


The  terms  on  the  right-hand  side  of  the  last  two  equations  of  (zo)  are  the 

rotatory  inertia  terms.  Since  the  beam  theory  assumptions  imply  that  the 

shear  resultants  Q  and  Q  are  zero,  the  usual  procedure  is  to  eliminate 
x  y 

them  from  the  first  two  of  (20)  by  the  last  two  equations.  The  result  is 


^y,  zz  Ry,  z  Px  pAu,  tt  P^yyU,  z2tt  +  *xyV>  **tt^ 


Mx,  zz  Rx,  z  Py  pAv,  p^xyu,  zztt  +  fxxvi  zztt} 


N  z  -  pz  =  pAw(  M 


(21) 


As  a  final  exercise,  the  moment  and  axial  force  resultants  are  eliminated 
from  Eqs.  (21)  by  means  of  Eqs.  (17),  with  the  result 

Ea  Ty,  zz  +  E(IyyU,  zzzz  +  V.  zzza>  +  pAu,  tt 


p^yyu,  zztt  +  *xyv,  zztt^  “  "Px 


(22) 


11 


o 


Ea  T  +  E(I u  +  I  v  )  +  pAv  „ 
x»  zz  xy  •  EZZS  XX  t  ZZZZ  t  tt 


^*xyu,  sett  +  *xxv,  sxtt^  "  "**y 

• 

t 

(23) 

-Eo  T  z  +  EAw(Zz-  pAw  tt  =  pz 

(24) 

( 

» 

t 

l 

f 


The  basic  equations  to  be  solved  are,  then,  Eqs.  (19),  (22),  (23),  and 

(24).  These  constitute  six  equations  in  the  six  unknowns  T,  T  ,  T  ,  u,  v,  w. 

x  y 

Note  that  only  the  first  of  Eqs.  (19)  and  Eq.  (24)  contain  w  and  T;  thus,  these 

two  equations  uncouple  from  the  rest.  When  I  =  0,  the  other  equations 

xy 

uncouple  similarly  into  pairs. 


J 
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SECTION  III 


ANALYSIS  OF  TIMOSHENKO  BEAM 

In  the  previous  section,  the  effects  of  rotatory  inertia  were  included 

9  10 

(hence  the  name  "Rayleigh  beam"  '  ).  In  this  section,  the  effects  of  shear 

deformation  will  be  included.  The  method  of  analysis  will  be  quite  analogous 
to  that  of  the  previous  section. 

The  basic  equations  (1)  through  (4)  obviously  will  be  unchanged.  The 
stress  resultant,  moment- resultant  equations  of  motion,  Eqs.  (6),  (7),  and 
(8),  will  also  be  unchanged.  The  heat  conduction  equations  (12)  and  (13)  will 
be  unchanged  except  in  the  term  involving  e. 

To  put  in  the  shear  deformation,  one  must  change  Eqs.  (14).  One 
assumes  that 


XI  =  u(z,  t)  ,  V  ~  v(z,  t) 

W  =  w(z,  t)  -  xpx( z,  t)  -  y0y(z,  t) 


(25) 


Here  0  and  0  represent  the  total  rotation  of  the  elements  including  the 
x  y 


rotation  due  to  bending  and  the  rotation  due  to  shear, 
imply  that 


The  assumptions  (25) 


'zz  =  E(w,  z  •  x  V  z  -  y*y,  z>-EaT 

e  =  (1  -  2v)(w  2-  x0x  z  -  y0y  z)  +  (3\  +  2p) 


(26) 


13 


If,  a*  before,  one  sets 


T  =  f  T  dA 
JA 


-J* 


dA 


T  =  f  xT  dA  ,  (27) 
y  •'A 


then  the  assumption  of  a  linear  temperature  profile  implies  that 


I  T  -  I  T 
T  T  +x  jcx_jr - xy_x 

*  I I -  I2 
xx  yy  xy 


I  T  -  I  T 

+  y  yy...x  -Jsy  j: 
ii  - 1 

xx  yy  xy 


(15) 


Using  (1  5)  and  (26)  in  the  heat  conduction  equations  yields 


k  T 


zz 


pc^l 


(3k  +  2^)Za2  Tq 


(k  +  n)pc 


2  2 


T  t  =  (3k  +  2p)a  T0(l  -  2v)Aw  zt 


kT  -  pc 
y,  zz  r 


1(  3  X.  +  2p)a 
1  (k  +  p)pc 


T  -k  A 

y.  t 


I  T  -  I  T 
xx  y  xy  x 

~I  I 

xx  yy  xy 


-(3k  +  2p)a  T0(l  -  2v)(Iyy^t  +  Ixy^it) 


kTx,zz-Pc  l  • 


(3k  +  2p)afc  Tf 


I  T  -  I  T 
k  A  -BL-* - 

II  -  I 
xx  yy  xy 


;  (k  +  p)pc  J  x,  t 

(3k  +  2p)a  Tq(1  -  2v)(Ixy^Xj  t  +  Ixx^y>  t) 


(28) 
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The  stress- resultant  deformation  equations  are  derived  as  before,  except 

that  the  Timoshenko  shear  constants  K  and  K  multiply  the  expressions  for 

x  y 

Q  and  Q  .  (For  a  more  complete  discussion  of  these  constants,  see  Ref.  11). 
x  y 

The  results  are 

N  =  EAw  -  Ea  T 
,  z 

M  =  -  E(I  e  +  I  e  )  -  Eq 
y  '  yyrx,  z  xyry,  z' 

Mx  =  ‘E(Ixy*x,  z  +  lxx+y,  z>  '  Ett 

Qx  =  VA(u,z-  *x> 

S  =  VA(V.  *  *  V 

These  expressions  must  now  be  substituted  into  the  equations  of  motion  (7); 
the  result  is 


(29) 


VA<V.»«-  V«’  '  py  *  pAv,  tt 

EA»,„-E«  T,,-  P,  = 

EV«.  -  +  Vy,  +  VA<“. .  ’  »«> 

-P'V„.«  +  Vy,«»+E“Ty,.  =  -Ry 

E^xypx,  zz  +  ^xx Py ,  ee'  +  EyPA^v,  s  “  Py^ 

-P<V*,t.+IxxV..'  +E«Tx.z  =  -Rx 


(30) 
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Equations  (28)  and  (30)  constitute  eight  equations  in  the  eight  unknowns  T, 

T  ,  T  i  u,  v,  w,  f  ,  f  .  The  boundary  conditions  are  the  usual  ones  for 
x  y  X  AT 

the  Timoshenko  beam1  z  plus  those  discussed  in  the  previous  section  for  the 
temperature  terms. 


SECTION  IV 


CALCULATIONS 


( 


n 


The  formula!  derived  herein  will  be  used  now  to  calculate  the  frequen¬ 
cies  of  the  free  oscillations  of  a  simply  supported  beam  for  which  I  =0. 

xy 

(This  is  not  necessarily  a  symmetric  cross  section.  )  Rotatory  inertia  and 
shear  deformation  will  be  neglected.  Two  cases  will  be  discussed:  axial 
vibrations  with  insulated  ends,  and  bending  vibrations  with  isothermal  ends. 
Bending  (and  temperature  variations)  will  be  assumed  to  occur  in  the  xz  plane 
only. 

With  the  above  assumptions,  the  basic  equations  become  - 


kT,zz‘pc 


l  (3K  +  2p)2a2  T0 

*  \  +  p. 


T,t  = 


(3\  +  2p)a  T0(l  -  2v)Aw  zt 


k  Ty,  zz  '  pc 


1  - 


(3\  +  2p)2a2  Tr 


TT 


T  -  2^  T 
^  Y 


-(3\  +  2p)a  T0(i-2v)Iyyu  zzt 


(31) 


El  u  +  pAu  +  Ea  T„  -  0 

yy  ,  zzzz  r  ,  tt  y,  zz 


EAw,  zz  -  pAw,  tt  -  Ea  T,  z  s  0 


(32) 


Thus  it  is  seen  that  these  equations  uncouple  into  one  pair  involving  w  and  T 
and  another  pair  involving  u  and  Ty.  The  equations  involving  w  and  T  are  the 
equations  of  axial  vibration,  and  the  equations  involving  u  and  Ty  are  the 
bending  equations;  they  will  be  treated  separately. 


17 


Case  I:  Axial  vibrations,  fixed  and  insulated  ends 


The  boundary  conditions  are 


w  =  T  =0  at  z  =  0,  ! 
.  z 


(33) 


The  basic  equations  are 


<3X  +  We  T„(l  -  2v  )Aw 
- ^ - --0 


w  — L  w  -  -T-  T  =0 

,  zz  2  ,  tt  A  ,  z 

C1 


(3X.  +  2(x)a  TQ 
*  (X.  +  p)pc 


2  _  E 

C1  "p 


(34) 


It  is  easily  seen  that  solutions  satisfying  the  boundary  condition  (33)  can  be 
taken  of  the  form 


w  =  F(t)  sin  ^ 

T  5  AG(t)  cos  — j— 


(35) 
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Then  Eqs.  (34)  become 


22  „„  <3h.  +  2p)a  Tft(l  -  2v)F 

(1  -  OG+2^-  G  +  ^ - - L  = 


2  2 

-L  F  +  2-V  F  -  G  =  0 
*1 


(36) 


Then  solutions  of  the  form  exp(iut)  exist  providing  u  satisfies  the  equation 


where  «  =  k/pc  is  the  thermal  diffusivity.  This  is  a  cubic  in  u.  If  <  =  0  or  is 
sufficiently  small,  there  will  be  two  roots  close  to  w  =  nircj/f,  which  corre¬ 
spond  to  the  free  vibrations  without  thermal  effects,  and  a  root  close  to 
«  =  i(n^w^/f  )[«t /(I  -  « )] ,  which  corresponds  to  a  decaying  "thermal  wave.  " 

It  is  of  interest  to  compare  the  oscillations  to  those  of  a  free  vibration  without 
thermal  effects.  One  sets 


w 


l 


and 


(38) 


Thus  (37)  becomes 


(;2  -  Hi  +  tvO  -  *&]  +  iV*  3 C  ■  0 


(39) 


while  <  is  usually  small  (<  a  10  ),  y  is  large,  at  least  for  the  smaller  modes, 

7 

(y  at  6  X  10  ).  With  this  in  mind,  one  finds  that  the  roots  close  to  (  *  1.  0  are 
given  by 


"  ‘  5^5tCK,2) 


h  V  rerfe*0*.  > 

*  l/\*  +  1  +  2,1 


The  root  near  £  =  i/y  is  given  by 


l  .i+ii-Uxi-L 
3  '  '  l/Al 


+  ex,2) 


Case  II:  Bending  vibration- -pinned  ends,  ends  kept  at  zero  temperature 
The  basic  equations  are,  from  (31)  and  (32), 


AStl-iT  A  T  ,  +  T0<*  -  2V>  -  u  J 

k  *  *Ty.  t  Ty  +  E  Vy  .  szt 


EI  u  +  pAu  „  +  Ea  T  _  =  0 

yy  ,  zzzz  r  ,  tt  y,  zz 


It  suffices  to  consider  solutions  of  the  form 


«/ .  t  »  ntr  2 
u  =  F(t)  Bin  — | — 


Ty  =  G(t)  sin 
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o 


Then 


(3X.  +  2|i)Q  Tq(1  -  2v)I 

K 


2  2 

yy  n  tt 


F  =  0 


G  =  0 


Solutions  of  the  form  exp(iwt)  exist  if 


As  before,  there  will  be  two  solutions  of  the  form 


(A)  = 


± 


(45) 


(46) 


and  one  of  the  form 


i 


(47) 


O 
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for  sufficiently  small  i .  We  are  mainly  interested  in  the  root  given  by 
Eq.  (46).  Accordingly,  one  sets 


O 
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be  it  temperature,  deformation,  or  street,  and  then  to  integrate  this  distri* 
bution  across  the  cross  section  to  obtain  resultants.  These  resultants  are 
then  the  quantities  which  are  considered  to  be  the  dependent  variables  of  the 
resulting  differential  equations.  Evidently,  the  success  of  the  method  depends 
upon  how  closely  the  assumed  distribution  approximates  the  actual  distribution. 
In  the  case  of  thin  beams,  plates,  and  shells,  a  linear  distribution  suffices, 
for  all  practical  cases,  to  describe  the  stresses  and  deformations.  For 
temperature  distributions,  the  situation  is  different.  Depending  on  the  prob¬ 
lem,  a  linear  distribution  might  or  might  not  be  indicated.  In  the  case  of  a 
plate  or  shell,  this  is  not  unreasonable;  in  the  case  of  a  beam,  the  success 
of  the  linear  approximation  definitely  would  depend  upon  the  individual  problem 
at  hand.  The  insulated  beam  was  chosen  for  two  reasons;  First,  it  is  believed 
that  the  linear  distribution  is  probably  a  good  approximation  here  since  there 
are  no  prescribed  temperature  gradients  over  the  cross  section.  Second,  it 
is  a  natural  boundary  condition  since  such  quantities  as  J 9T/8n  ds  vanish  for 
an  insulated  lateral  surface.  To  extend  beam  theory  to  include  the  other  types 
of  boundary  conditions,  one  must  use  some  further  approximations  to  obtain 
the  boundary  integrals  occurring  in  the  averaged  heat  conduction  equation. 

The  effects  of  shear  deformation  and  rotatory  inertia  are  definitely  of 

second  order  for  most  practical  cases,  becoming  of  importance  mainly  at 

high  frequencies,  or  where  the  shape  of  the  cross  section  changes  abruptly  or 

2 

rapidly.  The  Biot  coupling  is  also  a  small  effect,  becoming  of  importance 
mainly  when  the  dilatation  changes  rapidly.  This  would  more  or  less  imply 
that  for  dynamic  problems  both  effects  would  come  into  play  under  similar 
circumstances,  and  it  would  be  prudent  to  include  all  three  effects  in  any 
analysis  where  it  is  thought  that  any  one  effect  would  come  into  play.  (The 
case  of  a  short  beam  in  static  loading  where  shear  deformation  is  of  impor¬ 
tance  is  an  obvious  exception  to  the  above  statement. ) 
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